Sensitivity of solar f-mode travel times to internal flows 
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We compute f-mode travel-time sensitivity kernels for flows. Using a two-dimensional model, we show that 
it is important to account for several systematic effects, such as the foreshortening and the projection of the 
velocity vector onto the line of sight. Correcting for these effects is necessary before any data inversion is 
attempted away from the center of the solar disk. 
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I. INTRODUCTION 

The solar f-modes are useful tools for studying perturba- 
tions that reside within the top 2 Mm of the surface of the Sun, 
and have been used as a diagnostic as such before, e.g., Duvall 
& Gizon 1 2]. The linear forward problem of time-distance he- 
lioseismology is to calculate functions, or kernels, which give 
the sensitivity of travel-time measurements to these particular 
perturbations of the f-modes. For a discussion of this gen- 
eral procedure, which is the one that we will follow here, see 
Gizon & Birch [|. 

In this study we are interested inflow kernels, i.e., func- 
tions which give the two-dimensional sensitivity of travel- 
time measurements to small amplitude, two-dimensional, spa- 
tially varying flows. The flow kernel, K, is defined as 



Sr = JJ d 2 ru(r)-K(r), 



(1) 



where St is the perturbation to the travel time, u = (u x , u y ) is 
the two-dimensional flow perturbation, and the integral is over 
the position vector r = (x,y) over the whole solar surface. 
The kernel K = (K x , K y )isa two-dimensional vector which 
has sensitivity to flows in the x and y directions. 

We calculate the kernels from an approximate solar model, 
which utilizes the Born approximation in estimating the scat- 
tering from solar inhomogeneities. This has been shown to be 
more accurate than calculating kernels from a ray approxima- 
tion 1 1], since it includes finite- wavelength effects. 

Ultimately, the goal is to perform inversions with real time- 
distance measurements to extract the unknown u from Eq.{l}, 
given St and K. However, there are several difficulties in 
doing this since the travel-time measurements incur system- 
atic effects from the observing process, particularly due to the 
center-to-limb variations. To overcome these problems, one 
strategy is to calculate the sensitivity kernels from the outset 
accounting for the systematic effects so that the whole pro- 
cedure - the forward calculation and the inversions - remains 
consistent. 

We have accomplished this for two important effects: the 
foreshortening of the data pixel, and the projection of the so- 
lar velocity field onto the line of sight, both of which will be 
defined in detail below. We present an analysis of the contri- 
bution that these two effects have on the power spectrum of 
f-modes and on the kernels. For brevity, we leave the details 
of the full calculation of the power spectra and kernels for a 



future publication |4], and only focus here on the role of the 
systematic effects mentioned above on the final results. 



II. MODEL POWER SPECTRUM 

A. Line of sight projection 

A local region of the Sun may be approximated by a plane 
which is tangent to the sphere. The orientation of the normal 
to the plane with respect to the observer depends on the posi- 
tion on the sphere on which this plane lies. We denote the line 
of sight unit vector, i = (lh, i z ), as the direction which points 
towards the observer from the point where the plane intersects 
the sphere. In this notation, t z — 1 only at the disk center. 
Doppler measurements of oscillation velocities v on the Sun 
return the projection of v onto i. For the purpose of inver- 
sions, in order to have confidence in utilizing all of the data 
from instruments such as MDI and HMI, even data well away 
from disk center, the sensitivity kernels in principle should ac- 
count for a changing I, as to avoid the aforementioned inher- 
ent systematic errors. In the past, kernels have been calculated 
only with i z — 1 as a first attempt. We calculate kernels here 
for any i. 

In general, the observable oscillation signal, ip(r,t), can be 
written as 



ip(r,t) =&\l{r) -u(r,t)}. 



(2) 



The operator & involves the filtering procedure of the data, 
such as the time duration of the measurement, intrumental ef- 
fects such as the MTF and the foreshortening, and possibly 
some phase speed filter. Using a similar model for f-modes 
as in Gizon & Birch 1 3], and similar notation, we calculate an 
expression for v(r,t). The zero-order power spectrum of ip 
(in the absence of a flow) can then be written for horizontal 
wavenumber k and angular frequency uj as 



P(k,u) = (27r) 6 |£f(/c,cj)| 2 



(3) 



where 



uo) = (i£ h ■ k + y/m 8 (cj)G(k, cj)F(k, uo). 



(4) 
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FIG. 1: A cut in three different power spectra at constant uj/2-k = 3 mHz with B = 90°, where the x and y axes are in units of harmonic 
degree, (a) shows the power at the disk center = 0) for reference, (b) is the power at <I> = 75°, with the effect of the line of sight only 
(no foreshortening), (c) is again at <I> = 75°, but with the effect of foreshortening included (no line of sight effect). All figures have been 
normalized by their maximum values, and R s = 696Mra. We use a pixel size of an MDI hi-res image, dx = 0.43 Mm (see section llTBl . 



Here, k describes the frequency dependent damping, m s is the 
source power, G is the Green's function, and F is a function 
that describes the effect of the filter operator & mentioned 
above. For the purposes here, the important thing to note is the 
explicit dependence of the power on the line of sight vector, 
which is given as i = (— sin <1>, — cos 6 cos <£, sin 6 cos <£), 
where <£, the solar longitude, is measured from the central 
meridian, and 6 is the co-latitude measured from the solar 
north pole. In our local plane coordinates, is along x, and 
6 is along — y. 

Consider the power along the equator (<£ varies, 6 fixed at 
90°), and without any foreshortening effects for now. Upon 
close inspection of Eq.(|4}, one sees that the magnitude of the 
power along the k x direction is mostly unaffected, while the 
power along k y is reduced. This is shown in Fig.[T] where we 
plot cuts at constant angular frequency in the power spectrum. 
The first panel, Fig.[Tt, is at disk center, a reference plot where 
i — i z . Fig.[U> is a power cut at <£> = 75° towards the western 
limb, on the equator. As stated above, the power is reduced 
along the k y direction. 

To see this effect even more clearly, Fig.|2|is a plot of a cut 
in the power spectrum at constant uo for 4 different positions 
on the disk at fixed | k | , versus the direction of k (denoted by 
the angle k that k makes with the x axis), neglecting for the 
moment foreshortening effects. One sees that the power is 
always reduced along k y with respect to the value at the disk 
center. The power along k x is independent of when only 
line of sight effects are considered. 

Physically, we can understand this by considering the ex- 
treme case when we observe at the limb, on the equator. 
Waves that travel north- south will be undectable because the 
motion is perpendicular to the line of sight, while waves that 
travel along the equatorial plane still give power. 



B. Foreshortening 

Foreshortening occurs due to the fact that the spatial res- 
olution (measured on the Sun) decreases as we observe to- 
ward the limb. After projection onto the local plane, a ccd 
pixel will then appear elongated in the center-to-limb direc- 
tion. Again, consider a region towards the limb on the equator. 
A pixel that has a horizontal resolution of dx at disk center, 
now images a horizontal distance on the Sun of dx/ /i, where 
fji = cos sin O. This reduces the effective spatial resolu- 
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FIG. 2: Power in normalized units versus k for the 4 diffent values 
of $ along the equator (B = 90°), given in the legend. Here, £ = 
900 and lo/2tt = 3 mHz, corresponding to the power maximum 
of the f-mode. Note how only along k y (0 k = 90°, k = 270°) 
does the power get reduced due to the line of sight projection. Each 
curve is divided by the power at disk center (red dashed line) for 
normalization. 
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FIG. 3: Same as in Fig. with all curves at G = 90°, but now we 
show the power taking into account the foreshortening A(k) without 
any line of sight effects. Note how the foreshortening is 71/ 2 out of 
phase with the line of sight (compare with Fig.|2j. 

tion measured on the Sun in the center-to-limb direction, and 
hence leads to systematic observational errors. 

To correct for this, we return to the expressions for the 
power given above in Eqs. J3l4l) . The filter function F(k,u) 
contains the information about the foreshortening. We break 
the filter into its constituent parts, and rewrite it as 

F{k,u) = f{k,u)k{k), (5) 

where the lowercase f{k,uS) contains the MTF and/or other 
filters, which only depend on k, and A(fe) is the foreshort- 
ening correction. For the sake of simplicity, we obtain A(fe) 
by taking the Fourier transform of an elliptical disk, which 
is approximately the projection of a ccd pixel onto the lo- 
cal plane. We require that the ellipse only gets elongated 
along the center-to-limb direction, which is always along the 
semi-major axis. The semi-minor axis stays at fixed length. 
Carrying out this integration gives a term proportional to a 
Bessel function of first order. Expanding the Bessel function 
in a small parameter k (wavenumber) results in an expression 
given as 

A(fe) * Ao(fe) + Ai(fc) sin 2 (0 fe - a), (6) 

where Ao(fc) and Ai(fc) are only functions of k and are of 
little interest here, except to note that they depend on the ec- 
centricity of the ellipse, and hence on 6 and <£. We denote by 
a the angle which describes the orientation of the ellipse on 
the disk. More precisely, it is the angle that the semi-major 
axis of the ellipse makes with a local latitude line, and is also 
a function of 9 and All of the other details of this calcula- 
tion are beyond the scope of this paper, and will be presented 
elsewhere. 

The foreshortening term © now can be inserted into the 
power to study its effects. Returning to Fig.[T] in panel (c) we 
plot the power including the foreshortening function, but this 
time neglecting the line of sight effect. We see that the fore- 
shortening effect is 90 degrees out of phase with the line of 
sight effect from panel Fig.[U>. The power is reduced instead 



along the k x direction. We can understand this physically, be- 
cause waves are attenuated at high wavenumber, due to the 
smaller pixel size in Fourier space in the center-to-limb direc- 
tion. 

As was done for the line of sight case, we show a cut of the 
power 3t£ = 900, uj/2tt = 3 mHz versus 0k in Fig. [3]to see 
the reduced power in the k x direction more clearly as the limb 
is approached. This figure, along with Fig. |2] demonstrates 
that the net reduction of power is always larger due to line of 
sight effects, i.e., at $ = 75°, the power is only about 10 % 
of its disk center value along k y , while it is still about 70 % of 
its disk center value along k x . The evidence for this is also in 
Fig-ID^ where the power is almost gone along k y , but is still 
somewhat present in panel (c) along k x . 



III. KERNELS 

We calculate the flow kernels K based on the general 
method described in Gizon & Birch [3], where kernels were 
calculated only for damping and source perturbations to f- 
modes. The flow kernels here are computed in a completely 
analagous manner: we use the single- scattering Born approx- 
imation. The effects of the line of sight and foreshortening are 
explicitly included in the calculation. 

Since we are interested in detecting flows, we consider 
travel-time differences between two points on the Sun. Be- 
cause waves that travel in the direction of a flow travel faster 
than they would against it, this definition of travel time differ- 
ence is sensible for extracting flows. In Fig. 0] we plot kernels 
for travel-time differences that are sensitive to flows travelling 
in the +x direction (K x ), where the two observation points 
are given by the left and right crosses. Fig.|4t is computed for 
the center of the grid at disk center. For flow perturbations be- 
tween the observation points (blue region), one sees a negative 
value of the kernel, as expected. In other areas, wave effects 
produce non- trivial behavior. Fig.|4j), computed well towards 
the western limb, demonstrates how the line of sight and fore- 
shortening effects alter the kernel, as does Fig. 0J, which is 
a plot of a kernel on the central meridian at 6 = 15°. The 
kernel in Fig.0]; shows quite different structure, in fact, it has 
very little sensitivity along the line connecting the two obser- 
vation points; only waves that scatter off of the unperturbed 
ray path can be detected. The decomposition of the sensitiv- 
ity kernel into Fresnel zones (e.g. Jensen [5]) is not obvious 
in this last example. Also, note that even though the integral 
over space of the three kernels in Fig. |4]is approximately the 
same, the spatial-dependence is quite different, and an inver- 
sion using these kernels would give different results for flows 
varying on a scale comparable to or smaller than the wave- 
length of f-modes. 



IV. CONLUSIONS 

We have calculated sensitivity kernels for flows for time- 
distance helioseismology. The effects of the line of sight and 
the foreshortening, as defined above, have been demonstrated 
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FIG. 4: F-mode sensitivity kernels K x for flows in the +x direction, (a) is calculated at disk center = 0, G = 90°), (b) for $ = 75°, G = 
90°, and (c) for $ = 0, G = 15°. The crosses at n = (—5, 0) Mm, ri — (5, 0) Mm denote the two observation points. The units of the 
colorscale are s(km/s)~ 1 Mm~ 2 . Note that the integral over space of each kernel is approximately the same. 

for a model power spectrum and the kernels. The next step is spectra with the modeled ones to gain confidence in our theo- 
to study real data for f modes to compare the observed power retical sensitivity kernels. 
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